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I. INTRODUCTION 

Building analogies between fundamental field theories and condensed matter physics has been a long-standing trend. 
Usually, quantum field theory methods developed in fundamental interactions research were applied to condensed 
matter problems. However, a flow of ideas in the opposite direction exists as well. One of the fields influenced by this 
flow is the physics of strong interactions at low energies. As the perturbative QCD computation breaks down in this 
region, other methods are necessary. 

Savvydy [lj has demonstrated instability of the QCD vacuum with zero field strength in 1977. Since then exis¬ 
tence of various vacuum condensates in QCD became well established. This suggested looking for analogies with 
condensed matter models where such condensation occurs, superconductivity being the primary choice. Analogy with 
superconductivity gave rise to two models, one being the dual superconductivity 0 which is not our topic and the 
Nambu-Jona-Lasinio (NJL) model being the other. 

The Nambu-Jona-Lasinio model was first proposed in 1961 Q, well before QCD. Starting from the fermionic 
lagrangian built as an analogy to the BCS one 


C = iipdip + g 


(ipip ) 2 - (ip^y 


(i) 


Nambu and Jona-Lasinio demonstrated that fermion condensation occurs here as well. 

Since then, their model has been mostly used to describe mesons (see e.g. 0, 1). Attempts to look beyond the 
simple homogeneous case studied by the model’s authors were rare (one notable exception being !)• 

However, studying inhomogeneous cases is important as a number of interesting phenomena may occur. Consider 
the NJL gap equation with an Euclidean cutoff at p 2 = A: 


2t r 2 _ M 2 T A 2 ' 

^ _ 1 “ A 5 " n [F + 


( 2 ) 


This equation depends on M 2 only, which means it allows two symmetric opposite-sign solutions. 

This means that the model has two possible true vacuum states. Thus, there exists a possibility for existence of 
spatial regions with different vacuum values separated by domain walls. 

What happens to a fermion in such a vacuum? To answer this question, we note that as the condensate is obtained 
in the mean field approximation, the equation for fermionic states will be 


id — M (.x ) 


ip = 0. 


( 3 ) 


In the flat case where M depends on only one spatial coordinate this is exactly the Jackiw-Rebbi problem Q. In 
their 1976 work Jackiw and Rebbi have demonstrated existence of a E = 0 fermionic state with 1/2 fermion number. 
In our case where M(x) interpolates between this state’s asymptotic behaviour will be exp which 

means it will be bound to the domain wall. 

But what if the domain wall is not infinite but encloses a small bubble of —M vacuum in the +M sea instead? 
From general considerations we can expect that the bubble’s energy will be proportional to its surface area so it will 
not be stable but if there are fermions trapped at its surface their pressure may have a stabilising effect. If such stable 
bubbles exist, they are natural candidates for hadrons. 

Our paper is organised as follows: in section 2 we repeat our results first presented in Q where an approximation 
of the NJL gap equation that allows us to study inhomogeneous vacuum states was developed; in section 3 our 
solution to the Jackiw-Rebbi problem with a single flat domain wall solution of the approximate gap equation from 
[9] is presented and in section 4 we study the vacuum bubble and demonstrate that under certain conditions it gets 
stabilised. 


II. THE MEAN FIELD DYNAMICS 

The initial quark action of our model is 


S = 


ipidip + -^tp'y,_ l T a 'ipipj IJ -T a ip 


( 4 ) 
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Here ip bears the color index i as well as the spinor one and the generators T a of the SU ( N ) algebra are defined 
conventionally and obey the relations: 


rj-\arj~\b _ 


2 N, 


■S ab + - ( if abc T c + d abc T c ) ■ 


N 2 — 4 

^abc^abd c ^ad 

N r 


(5) 


N c is the number of colors. We keep it arbitrary for the sake of possible generalizations. 

The model 0 is not the Nambu-Jona-Lasinio model defined in j3]. However, its interaction is just the quark 
currents coupling and it’s intuitive to assume that if tha quarks are integrated out from the system somehow, then 
the resulting quark-quark interaction will look like this in the first order at least. 

We have proposed a derivation of such an action from the QCD lagrangian in 0, indeed. A similar quark 
interaction arises when instanton condensate is present (see El). 

The path integral of the model is equal to 


Z = / 'DQ'Dijj'D'ilj exp l—ild 


-Sp 


1 


7 „T a ( G12 + —ipip ) rT a ( G12 + —ipip 


ipidip + -^ip’y ll T a ipip^ ll T a ip- 


1 


2 G 


= J UQUipVip exp <^—i J d 4 x ipidip— 
—ip'y fl T a £lj fJ ’T a ip - G 2 5p7 M T a H7^T°H] } . 


( 6 ) 


up to a constant factor. The new variable H should have just enough degrees of freedom to couple to ipip in the action 
d6]). This means that H can be written as follows: 

n = w (a) T ( “ ) ; 

r («) £ |r( /3 ), t(^)t“| ; 
e {M7 5 ,y*,»y*7 5 } • 

The antisymmetric tensor section in 12 decouples from the theory because of the identity M = Au^v 11 

The Bogolyubov self-consistent field approach implies that we should study the lagrangian 


(7) 


C = 


ip (id - llM T a n c ^T a ) ip. 


( 8 ) 


where fl c is a solution of the classical equations of motion for the effective lagrangian obtained by integrating the 
quarks out of the action 0. The integration yields 


z = c vn 


exp {iG 2 J d 4 xSp-/ l ,T a njf 1 T a n} 


Det 


id — 7 M T a f 27 ^T a 


(9) 


which implies 


S e ff = iTr In 


id - - G 2 / d 4 x5'p7 M T a f2 7 /i T a f2. 


( 10 ) 


The corresponding equation of motion (the gap equation) is 


7 


id — 7 At T a f 27 At T a 


j^T c 


x ) - 2G 2 'y fl T a n(x)^T a = 0. 


( 11 ) 


If we assume 12 to be a constant scalar, then it’s just the case studied in [3], which corresponds to quark mass 
generation 7 Al T a 127 /i T a = M. But our goal is to study the inhomogeneous case. So we are going to put r ) fJ T a ^l^^T a = 
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M — Cl and expand the logarithm in (1101) around the constant solution. We obtain, omitting the constant part 
iTr Ini ^id — , 


= -G 2 J d 4 x&2 (0) + iSpS{ 0) J d 4 xtt(x) - 

~ ~ 2 S P J d 4 xd 4 yS(x - y)Cl(y)S(y - x)Cl(x) + ... ( 12 ) 

where S(x — y) is the plain Dirac propagator ( id x — M J S(x — y) = S(x — y) and we have defined 


Sp^^T a n{x)^T a VL{x) = $ 2 (fi). 


(13) 


The first term in the expansion can be rewritten as follows: 


s pS(0) j d 4 xti(x) =Sp j d (2*)? (M) 


Let us now demonstrate that each term in the series can be reduced to 


d 4 pd 4 x^Ej(x,p). 

i=0 


For this, it’s enough to demonstrate that 

J d 4 yS(x — y)£l{y) J d 4 pe~ ip( - y - z) E{p, z) = J d 4 p e - ip(x " z) S(p, z). 


(15) 
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So 


d 4 yS(x — y)£l(y) J d 4 pe ip( - v z) E(p,z) = 
f d 4 yd 4 gd 4 p __ ia(x _ v) q + M - 


(2tt)4 

.;4„,j4„j4„ 


q 2 - M 2 
f d 4 yd 4 qd 4 p _ iq(x _ y) q + M 

J (2tt) 4 q 2 - Af 2 


n(i/)e- < p("-*)s( P ,«) = 


E 

71—0 


(y-z) ai ~(y-z)° 


n\ 


d d 


dz a i "dz“’ 


-S2(z) 


3 -<p(l/-*)g(p j *) = 


[ JPy&qtPp - iQ(x - v) q + M 

J (2t r) 4 q 2 - M 2 


E; 

71=0 


9 a 


H) 


n! [az" 1 az“" j L^Pai ap, 


a a 


0 -ip(y~ z ) 


l(p, Z) = 


d yd qd p lq(x _ y) q + M _ ip(y _ z) 
(2t r) 4 q 2 -M 2 


(~i) r 


n\ 


a 


a«“! ”az“ 


5 n(z) 


><E 

71=0 

= [ d 4 qd 4 pe-^~^-^- 


a a 


apai "d Pl 
f d 4 y ; (v _ 


(P, z) = 


E 

71=0 


H) T 


q 2 -M 2 7 (2tt) 4 

a a 


e i(p-q)y x 


a 3 


! La^ ai ''az a " v 'J dp ai "dp, 


i(p, z) = 


= J dVpe-^-^l ±iL 


S(p-q)Y 


nr 


E 


nr 


71=0 

' a 


3 3 -nw 


! La^ ai "az“» J dp ai "dp, 


a 


a a 


0, z) = 


= f d 4 vc- ip(x ~ z) P + M 

J p 2 — m 2 ^ n! I az ai az“ n v ~'j a PQ1 ap Q 

n=0 - J 


a - 

-ft(z) 


a a 


HH> Q.E.D. 


Therefore, it’s now obvious that if we define the operator 


V(x,p) = Y 


nr 


71=0 


a 


5 -f 2(z) 


! [dz ai "dz a " v ’\ dp ai "dp a J 


a a 


it’s possible to express the expansion m as 


S eff [fi] = -G 2 / (ft) + i 


d 4 xd 4 p 


^E 


(-i) 


fc-i 


fc=i 


p + M 
k [p 2 - m 2 




o 1. 


We can easily read from this that all the integrations in (fl2l) reduce to expressions of type 


I(m,n) ai .. am = J d 


4 Pctl”PcXrt 

Pi 


{p 2 - Af 2 ) ' 

First of all, it’s easy to see that integrals like this are nonzero for even values of to only and are finite if 

d(n, to) = 2n — to — 4 > 0. 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


It’s obvious that differentiation of the integrand in (fl9l) increments its d(n, to) by 1 while leaving the difference 
n — to unchanged. Thus we can conclude from (1181) that the series contains only the integrals (1191) with m < n. But 
it’s also true that a fc-th order in Q term can contain only the integrals with n > k. 

Therefore we immediately conclude that the expansion (1181) has only finite number of divergent terms which are 


1. first order term; 
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2. second order up to second derivatives of Cl; 

3. third order up to first derivatives; 

4. fourth order with no derivatives of Cl. 


Regularization by cutoff will replace the divergencies with some finite factors of order A n in the cutoff parameter. 
Therefore, the finite part can be treated as a small correction to the divergent one. So by omitting the finite terms we 
can build an approximation of the effective action dll that has remarkable features. First, it’s similar in structure 
to the Ginzburg-Landau functional in superconductivity, since it contains derivatives up to second and nonlinearities 
up to fourth order; second, it requires no extra conditions being imposed on the model to be valid. 

To calculate the terms of our interest we use the Passarino-Veltman reduction fl2| 


= 0 ; 


/(2m T 1, n) CK1 .. ( x 2m _j_ 1 

/(2rn, n) ai ..a 2m = CS (?7a 1 a 2 ■■Va 2 m-i 0 ‘ 2 m) i J d p 


[p 2 - M 2 ] 7 


( 21 ) 


where C is a symmetry factor and S (v ai a 2 -- r lot 2 m-iot 2 m) is a symmetric tensor power of the Lorentz metrics. 

After rather long computation we obtain the following expression for the approximate ’’Ginzburg-Landau” action: 


= -G 2 f d 4 x<S> 2 (Cl) + 

(K + l) 

' f d 4 xSp {2 M [Z( A) - A/ 2 ] Cl (x) + M 2 [n(®)] 1 


In 


Z( A) - 2 M 2 


j l _ l Cl(x) r y ,J 'Cl(x) + MCl(x)idCl(x)— 


-i dCl(x) --^p ll (d v Cl(x))^d l 'Cl(x)- 
-^-Cl(x) lll Cl(x)-i ll Cl(x)+ 

+ ^ 1 ti Cl(x)Y [d„f2(a;)] YCl(x) ~ ^7 M H(a:)7 M n(a’)3f2(a;)-t- 

+ 777 Sp [y l j.Cl(x)'y , 1 Cl(x )] 2 + ■^-Sp"f ll Cl(x)'y l 2Cl(x)'y IJ ’Cl(x)'y''Cl(x) 
z4 4o 


( 22 ) 


where it has been put 


Z(A) 


A 2 


(23) 


Note that differential and nonlinear terms are of the same logarithmic order in A. This implies that nonlinearities 
play an important role in the inhomogeneous case and can’t be treated as just small corrections. 

Let us now change variables back by substituting Cl M — Cl. In the fourth order we get 


Sp"/ /J .Cl"f^fl'y ll fl^ v n Sp 

-SMCl^nrCl + 4:M r j ll Cl r y IJ, Cl'y v Cl'y 71 — 8AI 2r y ti Cl'y^Cl+ 

+16 M 2 Cl 2 - 2 M 2 ( 7 ^ 7 ^) 2 + 32 M 3 Cl - 8A/ 4 ] ; 

Sp [y l j,Cl(x)'y IJ Cl(x)] 2 Sp ['y^Cl(x)'y> J -Cl(x)] 2 + 

+Sp [—8A iCl^^Cl^Cl — 2M r y IJl Cl r y fJ 'Cl^ 1 2Cl r y v + 

+16M 2 Cl 2 + M 2 (j^Clj 17 ) 2 + SM^Cl^Cl- 

-64 M 3 Cl + 16A/ 4 ] . (24) 


The third order yields 
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SpCl'y fl Cl / y^fl —> —Sp 

- 4 Mn 2 +12 M 2 n - 4M 3 ] ; 

Sp 7 M fi(a;) 7 M [<9„fi(a:)] 7 "Cl{x) — ^ /Ji Q(x)'y IJ 'Cl(x)dCl(x) —» 

—Sp 7 M f 2 (a;) 7 ' i [<9„fi(a:)] ^Cl(x) — 'y /J ,£l(x)'y IJ ‘£l(x)d£l(x)+ 

+M'y li Q'f [ 7 ", &,n] + 4 Mfi [ 7 ", &,n] - 4 M 2 [ 7 y , a,n]]. ( 25 ) 

Finally 

Spfi 2 -> Sp [fi 2 - 2MQ + M 2 ] ; 

Spy^fiy^fi —> Sp [y^Uy^fi — 8MCl + 4M 2 ] . (26) 

Therefore, after omitting the constants and subtracting full derivatives where appropriate, we obtain 

S = —G 2 J d 4 x $ 2 (fi) + 

In (^ 2 + 1 ) /■ r zfA) - - M- 

H-—- J d 4 xSp^ -^7 /x n(a:)7 /i n(a:) - y fl(a:)iaf2(a:)- 

—g-7^7 # ‘ - - [an(a:)J - y 7 M (0„fi(aO) 7^0(3;)- 

-y M fi(a;)7 A1 [c^f^a;)] Y&{x) + y M f2(a:)7 M f2(a;)afi(a;)+ 

T-^p [7^0*07^ 0*0] 2 + ^Spy M f2(a07^(z)7^(z)7^(z) j • (27) 

Now we note that 

fi = 0 + 7 m w^; 

0,W„€ { / + F a T a + iy 5 (5 + G a T°)| f,F a ,g,G a £ R} . (28) 

Then it immediately follows from the 7 -matrices’ properties that 

7 M 0 = </) + 7 /i ; 

7 M wy = w+ 7 y 

7^7^ = 4</> + - 2 7 ^W+. (29) 

Using the identities above we can easily demonstrate that 

f d 4 xS P ^-n(x)idn(x) - [-/\dM] =0 (30) 

and the action becomes 

S = -G 2 J d 4 x<S> 2 (fi) + 

+ ln ^: - / ^p{-^7^)7^)- 

-i pU(x)] 2 - y, (^fi(*)) 7^U(*)~ 

-y M U(;r) 7 M [9„n(a:)] y^a;) + y /i U(a;)7' i O(a;)9U(a:)+ 

T^Sp [7M^(*)7 #i ^(*)] 2 + ySpy M f2(a07^(z)7^(z)7^(z) j . (31) 



Despite the simplification, the dynamics defined by m is still very complicated due to the big number of degrees 
of freedom. That’s why we are going investigate simple special cases. First of all, it’s obvious that we can drop the 
color sector putting </> = £ + + i^w^. The next question is whether it is possible to put W M = 0. 

The only nonvanishing terms in (1311) are those that contain even numbers of 7 -matrices. Thus, the only terms linear 
in Wp that would produce a nonvanishing part of the equations on come from the third-order first-derivative part. 
However, it’s easy to demonstrate that these terms are equal to 2 iSp [{W^ + — bF M ) and the corresponding 

part of the motion equation vanishes if (f>(x ) = £(x) e R. 

The action for £ is 


+ 


= —G 2 J d A x^ 2 (x) + 
(& + i) 


N r . In 


8tt 2 


d A x { -Z{A)e(x) - + ^ 4 ) \ = 


= -K J d A x | \d^d^ + 


Z{ A) 


87 r 2 G 2 


,2 c 4 


r- 


(32) 


^ln(^ + l)_ 

or, if we take into account that M is obtained as a constant solution of the gap equation ED, we can finally see that 

16t r 2 G 2 


: [£] = K f d A x | + 


M 


N c In 


A 2 

TP 


1) 


e-t ■ 


(33) 


III. THE BOUND STATES AT A DOMAIN WALL 


The (1551) is nothing other than the real </> 4 model. The equation of motion is 

□£ - 2A 2 £ + 2£ 3 = 0, (34) 

where we have denoted A 2 = M 2 -|- 16 7 —r. It has the well-known kink solution 

£(x) = Atanh [Ax 3 ] . (35) 

This solution interpolates between £ = —A and £ = A. These are the two vacuum states of the action (1331) . Thus, 
our solution describes a domain wall between regions of space with different vacua. 

Let us now solve the Dirac equation with (j> A kink potential. The ’’wideness” parameter of a kink can be absorbed 
into its ’’height” by scaling the spacetime variables appropriately: 

r = Ax°; 
z = Ax 3 . 

(36) 

Thus the equation becomes 

[*7°i9 r + *7 3 i9 2 - /itanh ( 2 )] ^(r, z) = 0, (37) 

where /i = 1 in our case but we are going to keep /i as an arbitrary parameter for a while. 

The next step is putting 


= e~ iEr 


(f>(z) 

X(z) 


(38) 


By choosing the appropriate Dirac matrix representation we reduce the equation to the following system: 

ifJz-r 1 + [E — /utanh(z)] <p = 0; 
dz 

ivz^r- + [E + /rtanh(z)] % = 0. 
az 


( 39 ) 
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Then we make an ansatz 


cj>(z) = <j>{z) |t); 

X(*) = x(*)|t>, (40) 

where u z |f) = |t) ( we might have chosen the other eigenvector of the Pauli matrix which wouldn’t change the picture 
much). 

The resulting system of equations for <fi, \ can be written in the following matrix form: 


d (f> 
dz X 


[iEcr x — /itanh(z)(jj / ] 


X 


(41) 


Now let us define 


These vectors satisfy the identities: 


We can now make the following ansatz: 

.X. 

which leads us to 


1 + ) 

b> 


1 

71 

i 

71 



a v 1 +) — 1 +); 

°y b> = - |->; 

Vx 1+) = i b); 
°x b) = -t|+) 


= p(z) cosh M ( z ) | + ) + lb) cosh* 1 ( 2 ) |—), 


(42) 


(43) 

(44) 


d 

p{x) 


dz 

q{x) 



0 Ecosh 2fl {z) 

-Ecosh.~ 2fl (z) 0 

Then we exclude q(z) from this system of equations and obtain the following: 

cosh _2Al (2) dp 


p(x) 

q{x) 


(45) 


q(z) = 


E 


dz 1 


- 2/xtanh (z)^p + E 2 p(z) = 0. 


(46) 


If E = 0 p and q decouple from each other and we can just put p = 1, q = 0; the other linearly independent solution 
diverges at the infinities and we don’t take it into account. 

Finally, we put £ = sinh(z) and the last equation turns into 


(1+< j )|| + (1-2 mk| + E 2 p(O=0. 


(47) 


This equation falls into the hypergeometric class, furthermore, it can be easily proven that there’s a series of polynomial 
solutions that are orthogonal at (— 00 , 00 ) with the measure 


Their spectrum is 


^(C) = (1 + C 2 ) 11 ^ 

E 2 n =n (2/z - n). 


(48) 

(49) 


It can be obtained by studying the asymptotics of the equation and taking into account that the leading order should 
disappear. These solutions with \n\ < p correspond to the bound states. The Rodrigues formula for these polynomials 
is 


i>„«) = (i+<y +4 |b+'; 2 r 


(50) 


One can note that the polynomials are very similar to the well-known Gegenbauer polynomials. See the_appendix for 
more details. For the complete analysis of the equation (1471) in terms of hypergeometric functions see 13]. 
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IV. SOAP BUBBLE HADRONS OR BOGOLYUBOV BAG? 


The equation (l34l) has two possible vacuum states ±A. The kink solution (l35l) describes a domain wall between 
the different-vacuum regions of space. It depends on one parameter A which defines both vacuum condensate and 
thickness of the wall, the latter being proportional to A 

Let us now consider a case when a spherical bubble of radius R of—A vacuum is trapped inside the +A-vacuum 
space region. If the bubble is large enough (i.e. R A -1 ), the domain wall at its edge will behave almost as in the 
flat case studied above. This means that there must exist a quark state bound to the bubble surface. 

The stress-energy tensor corresponding to the action (1331) is given by the formula 


r p^ v = I\ 




d^d v i - - rr a 2 £ 2 + -ve 


VT ,a 


(51) 


Energy density is given by its component T 00 . For a stationary configuration, taking into account the equation of 
motion we obtain simply 


£ = 




(52) 


If we have just a locally-inhomogeneous solution that approaches vacuum value at infinity, such as our domain wall, 
the zero energy level must correspond to vacuum, so the energy will be given by 


E = ^Jd 3 x [A 4 — £ 4 (a’)] . 

Substituting the kink solution into it and integrating over the ^-coordinate yields 


(53) 


E = J KA 3 dS, (54) 

where the integration is carried out over the domain wall surface. 

This means that the energy of a large spherical vacuum bubble will be AnKA 3 R 2 and the bubble must be unstable. 
Such a bubble would collapse spontaneously. 

However, if there are quarks bound to the bubble surface, the picture changes drastically. When the bubble’s radius 
becomes comparable to A^ 1 , its boundary can no longer be considered thin; the bubble’s energy gets distributed over 
the whole its volume and the bubble no longer contains —A-vacuum. Instead, it can be approximated by a bubble of 
£ = 0 inside the £ = A vacuum. 

Such a bubble containing quarks is exactly the Bogolyubov bag (see 0 for details). The quarks in such a bag have 
energy proportional to R~ 3 . So, the total energy of the system will have the form aR 3 + ^ and have a minimum at 
some nonzero value of R. 

An interesting question arises here: whether a bubble indeed collapses into a Bogolyubov bag or it is stabilised by 
quark pressure at some equilibrium radius. 

First of all, it is obvious that to answer this question, we must study a bubble with R A -1 . Thus, we can treat 
the boundary as almost flat and solve the problem perturbatively. 

Let us write the Dirac equation in spherical coordinates following the notation of 0). The Dirac spinors will have 
the form 




K, 


9{ t )Xk 

f{r)x 


(55) 


where are eigenvectors of the operators K = /3 ycrl + lj and z-component of the full angular momentum j = 
Then the Dirac equation will yield 


[E - M{r))g = i^ + 

[E + M(r)} f = ij- + 

dr r 


r 


( 56 ) 
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Now we can define 


which leads us to 


g{r) 

fir ) 



Eu = M(r)cr z u 


du K 

+ i<Jx~r + -v v u. 
dr r 


(57) 


(58) 


As R A -1 , we can approximate M(r) as follows: 

M ( r ) = A tanh [A (r — i?)]. (59) 

We are interested in the behaviour near r = R. Thus we can assume 1/r = 1/R and treat it as a constant in (15811 . 
Making the same ansatz as in (l44l) then leads us to the equation 

||-2ta„hw| + (jS 2 -^)p( 2 ) = 0, (60) 

which differs from (1461) only by an addition to E 2 which means simply a shift of energy levels. For the zero mode we 
obtain 


E = ±^. (61) 

The proper perturbative treatment of the problem, however, requires us to take into account corrections to (1591) . 
We can build the solution to the radial version of the <jA equation (1341) as a power series in 1/R with (15^1) as the zero 
order. Obviously, this series will have a nonzero first order term which must be taken into account to calculate the 
exact 1/R correction to the zero mode energy. 

We will not go into the details of exact computation. It is enough for us to know that the energy of the bubble 
with a trapped quark is 


E = 4nKA 3 R 2 + —, (62) 

R 

where v is a constant not depending on R. The energy (1621) will have a minimum at 



This radius depends on the constant K which, in turn, is implicitly related to the NJL mass gap M. Let us study 
this relation in more detail. 

We regularize the equation m by an Euclidean cutoff at p 2 = A 2 to arrive at 


A 2 - M 2 In 



8t r 2 G 2 
N c 


(64) 


As the fourth-order term in ([4j) is proportional to G -2 , we will study the small G case G< A which corresponds to 
strong interaction. After expanding the logarithm in (1641) into a power series we see that in the small G case -A- must 
be small as well. Thus, for an approximate value of M we can leave only the first nonvanishing term in the expansion 
which yields 


v^A 2 




47rG 


(65) 


This, together with (l33l) . leads us to the conclusion that with the same degree of accuracy K = This means, 
that, if the interaction is strong enough, we can expect the equlibrium configuration of a system of quarks trapped in 
a vacuum inhomogeneity to be a ’’soap bubble” - in the sense that its energy and particle density is concentrated at 
the surface. 
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V. CONCLUSIONS AND DISCUSSION 


We have developed a Ginzburg-Landau-like approximation of the NJL gap equation that allows us to explore 
inhomogeneous vacuum configurations. We have demonstarted that this approximate equation of motion has a simple 
scalar sector which is nothing other than the (p A model. It is well-known that the phi 4 field equation has the kink 
solution which in our case has the physical meaning of a domain wall between space regions with opposite quark 
condensate values. 

In such an inhomogeneous vacuum Dirac equation for quark states can be solved exactly. The general problem 


d — /xtanh(z) 


ip = 0 


( 66 ) 


has a series of bound states with energies given by 

En =n (2/r — n ); \n\ < //, (67) 

however our case is equivalent to p = 1 which means that the only bound state has E = 0. 

if the domain wall is not infinite and planar but encloses a ball of opposite-sign vacuum, this bound state energy 
gets an 0(R~ 1 ) correction which, in turn, implies that under certain conditions the system can be stabilised by quark 
pressure resulting in a peculiar “soap bubble” configuration which has its energy and matter concentrated at the 
surface. 

Such “soap bubble” states are natural candidate for hadrons. However, such a model for hadrons fails to explain 
the hadron mass hierarchy. 

If quark pressure and surface tension are the only forces thet define bubble size, then the more quarks are trapped 
in the bubble, the bigger and, as the domain wall energy is proportional to surface area, heavier it must be. 

Thus, the model predicts existence of a single-quark configuration lighter than a proton which is not observed. 

So, modifications to our model are required to explain the real hadrons. One possibility for such a modification is 
some residual QCD attraction between trapped quarks. If such an attraction is strong enough, then, as the quarks 
are light, the single-quark configuration will become unstable due to pair production. 

Different bubble surface geometry is another possibility. Strictly speaking, sphericity of the bubble in case of several 
quarks present has yet to be proven. It is worth noting here, that if the two-quark bubble is indeed not a sphere 
but an elongated ellipsoid stretched along the line that joins the quarks, it will produce the effective linear potential 
which is expected for confinement. 


VI. APPENDIX. THE FERMIONIC STATES 

The bound states’ wavefunctions should be orthogonal: 


J d 3 xip+ ( x , y, z)ip m {x , y, z) = N(n)S nm . 


( 68 ) 


In our case the states are actually ’’semibound” since the fermion motion is restricted along the z-axis alone. Thus 
the normalization constant will be infinite. So we actually should prove that this scalar product reduces to 


J d 3 x*P+(x,y,z)iPm(x,y,z) = J dxdy J dzip+{z)ip m {z) 


(69) 


and the z-integral is actually finite and J dzip+(z)ip m (z) = C n 5 nm for the wavefunctions that don’t depend on x,y 
that were studied above. 

To achieve this we should first reconstruct the complete Dirac spinors for the (l50l) solutions. Taking into account 
the transformations from section 3 we obtain 


Vvr ( z ) 


cosh M (z) 

V2 


(p„(sinh(^)) + J- ^^ h(z)) ))lt) 

t(p n (sinh(«))-^^^i))|t) 


(70) 


Thus the normalization condition becomes 
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(nt|wt}= J dxdy J dz cosh 2fl (z) x 

1 dpn(sinh(z)) dp m (smh(z)) 


Pn(smh(z))p m (sinh(z)) 4 

E n E m dz dz 

However, for the bound states we have, by substituting again sinh(^) = £ 


(71) 


dzcosh 2fJ, (z) 


Pn (sinh(z) )p m (sinh(z)) + 


1 dp n (sinh(z)) dp m (smh(z)) 


E n Er, 


dz 


dz 


+oo 


= J dc(i + cT^ 

—oo 

+oo 

= [ dca+cr^ 


Pn(C)Pm( C) 


1 + C 2 dp n {£) dp m (() 


E n E m d( d( 


, n 1 - dpn(C) 1 + C 2 d 2 p n { C) 

Pn[U E n E m ^ d( E n E m dC 2 


Pm{ C) = 


En 


+oo 


= [l + ^) j d(( 1 + C 2 ) 2 PniOPmiC)- 


(72) 


Therefore we can conclude that solutions with E n = —E m (we can make a convention E_ n = — £ , „,p_„(£) = p n (£)) 

+oo _ _ 1 

are orthogonal. For the rest of the cases scalar product reduces to / d£ (1 + £ 2 ) M 2 Pn{()Pm(C) an d we can see 

— CO 

that the measure is exactly (H51) . 

Only the opposite spin case remains now uninvestigated. However orthogonality of the opposite spin solutions is 
guaranteed trivially by (f | 4-) = 0 and the spin-down solution are very similar to the spin-up case studied above, 
that’s why we aren’t going into studying them in detail. 
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